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1 Introduction 

In the present article we will focus on developments concerning harmonic and 
other moments of a domain in the complex plane, the exponential transform 
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and, in the case of a quadrature domain, the relation between these objects 
and the rcsuhant and ehmination function. Much of the material is based 
on joint work with Vladimir Tkachev, Mihai Putinar, Ahmed Sebbar and is 
in addition inspired by ideas of Mark Mineev-Weinstein, Paul Wiegmann, 
Anton Zabrodin and others. 

The organization of the paper is as follows. Section 2 gives the basic defi- 
nitions of Cauchy and exponential transforms, including extended versions in 
four complex variables. In Section 3 it is shown how these generate harmonic 
and exponential moments. In the case of quadrature domains (equivalently, 
algebraic domains or finitely determined domains) , there are strong algebraic 
relationships between the various transforms and moments, and this is ex- 
posed in Section 4. Most of the material in Sections 2-4 is by now relatively 
classical, being based on developments in the period 1970-2000 within com- 
plex analysis and operator theory. For example, the exponential transform 
came out as a byproduct from the theory of hyponormal operators. 

Central for the paper is Section 5 on the meromorphic resultant, which 
was introduced in [19]. One main message is that, in the case of quadra- 
ture domains, the elimination function, which is defined by means of the 
meromorphic resultant (and hence is a purely algebraic object), turns out to 
be the same as the exponential transform (which is an analytic object). In 
Section 6, some potential theoretic interpretations of the resultant are given. 

In the last three sections. Section 7-9, we connect the previously dis- 
cussed material to some quite exciting and relatively recent developments 
in mathematical physics. It much concerns deformations (or variations) of 
domains when the harmonic moments are used as coordinates (in the sim- 
ply connected case) . Such deformations fit into the framework of Laplacian 
growth processes, or moving boundaries for Hele-Shaw flows. More exactly, 
they can naturally be thought of as evolutions driven by harmonic gradients. 
We mostly restrict to finite dimensional subclasses of domains, namely poly- 
nomial (of a fixed degree) conformal images of the unit disk, this in order to 
keep the presentation transparent and rigorous. A nice feature then is that 
the Jacobian determinant for the transition between harmonic moment co- 
ordinates and the coordinates provided by the coefficients of the polynomial 
mapping functions can be made fully explicit, (7.5), in terms of a resultant 
involving the mapping function. This result is due to O. Kuznetsova and 
V. Tkachev [25], [46]. 

The further results concern the string equation (7.8) (which is equiva- 
lent to the Polubarinova-Galin equation for a Hele-shaw blob), integrability 
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properties of the harmonic moments (3.6), and a corresponding prepotential 
(3.8), which is the logarithm of a r-function and which in some sense ex- 
plains the mentioned integrability properties. Finally, we give Hamiltonian 
formulations of the evolution equation for the conformal map under variation 
of the moments (9.8). All these matters are due to M. Mineev-Weinstein, 
P. Wiegmann, A. Zabrodin and others, see for example [53], [29], [27], [24], 
but we try to explain everything in our own words and using the particular 
settings of the present paper. 

1.1 List of notations 

Below is a list of general notations that will be used. 

. D = {C e C : Id < 1}, D(a,r) = {( e C : \C - a\ < r}. 

• P = CU{oo}, the Riemann sphere. 

• dm — dm{z) — dx Ady — ^^dz A dz. {z — x + iy). 

• /*(c) = /(i/d- 



2 The Cauchy and exponential transforms 

The logarithmic potential of a measure /i with compact support in C is 

U^{z)^ j log^l^d//(C), (2.1) 



satisfying — Af/^ = 27r/i in the sense of distributions. The gradient of f/^ can, 
modulo a constant factor and a complex conjugation, be identified with the 
Cauchy transform 

2dU,{z) 1 fdfiiO 



n dz n J 



which may be naturally viewed as a differential: C^{z)dz. Prom these po- 
tentials, or fields, /i can be recovered by 

^ dz 2Tr ^ ^ ^ 
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We shall mainly deal with measures of the form /i — Xn?Ti, where m 
denotes Lebesgue measure in C and O C C is a bounded open set. Then we 
write [/n, Co in place of U^, and C^. On writing Co as 



C^iz) = — f 
^ ^ 27ri 7n 



- AdC. 



C-z 27ri C 



one is naturally led to consider the more symmetric "double Cauchy trans- 
form" , 

Ca{z,w) = ^ [ -^A J^, (2.3) 

which is much richer than the original transform. It may be viewed as a 
double differential: C^{z,w) dzdw. After exponentiation it gives the by now 
quite well studied [5], [34], [36], [37], [13] exponential transform of ft: 

Ea{z,w) ^ expCn{z,w). 

The original Cauchy transform reappears when specializing one variable at 
infinity: 

Cci{z) — Res Cci{z,w) dw — — lim wCq{z,w). 

■w=oo w—>-oo 

As a substitute for (2.2) we have, for the double transform, 

d^Cn{z,w) 



dzdw 



-nSiz - w)xn{z)xn{w). 



Even the double Cauchy transform is not fully complete. It contains 
the Cauchy kernel which is a meromorphic differential on the Riemann 
sphere with a pole at ( — z, but it has also a pole ai ( — oo. It is natural to 
make the latter pole visible and movable. This has the additional advantage 
that one can avoid the two Cauchy kernels which appear in the definitions of 
the double Cauchy transform, and the exponential transform, to have coin- 
ciding poles (namely at infinity). Thus one arrives naturally at the extended 
(or four variable) Cauchy and exponential transforms: 

Ca(.,w-,a,i) = ± [(^-^)A(jL.JLv (2.4) 
STTiinC-^ C-« C-w C-b 

En{z, w; a, b) = exp Cn{z, w; a, 6) = , , , r. (2.5) 

EQ[z,b)EQ{a,w) 
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Clearly Cq{z,w) = Cq{z,w;oo,oo), Eq{z,w) = E fi{z, w; oo, oo). If the 
points z, w, a, b arc taken to be all distinct, then both transforms are well 
defined and finite for any open set Q in the Riemann sphere P. For example, 
Ep{z, w;a,b) turns out to be the modulus squared of the cross-ratio: 

Ev{z,w;a,b)^\{z:a:w:b)\'^, (2.6) 

where, according to the classical definition [2] , 

(z — w)(a — b) , ^ 

{z:a:w:b) = ) ( 2.7 

[z — b){a — w) 

with the variables in this order. To prove (2.6), one may use the formula for 
the two- variable exponential transform for the disk Q = ©(0, it!) when both 
variables are inside, namely (see [13]) 

2 



^D(o,R)(^,w) = ' _ - {z,w eB{0,R)), 



insert this into the last member of (2.5) and let R — >■ oo. A different proof 
of (2.6) will be indicated in Section 6. 

We record also the formula for the exponential transform for an arbitrary 
disk D(a, R) when both variables are outside {z,w e C \ D(a, R)): 

^B(a,.)(.,-) = l-^^3^^--y. (2.8) 

See again [13]. 

From (2.6) it is immediate that for any domain ^2 C P, 

EQ{z,w;a,b)Ep\a{z,w;a,b) = \{z : a : w : b)\^ (2.9) 

(clearly Ed makes sense even if the set D is not open). The two- variable 
exponential transform Ep\q{z,w) is identically zero if fl is bounded, but 
there is still a counterpart of (2.9) in this case. It is 

En{z, w) ■ / ^\z- (2.10) 

where Hq{z, w) is the interior exponential transform^ which is a renormalized 
version of one over the exponential transform of the exterior domain, more 
precisely 

Hn{z,w) = lim -— . 

R-s-oo -rt^-C/D(0,il)\n 
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We consider this function only for z,w E Q, and then it is analytic in z and 
antianalytic in w. 

When the integral in the definition of the exponential transform is trans- 
formed into a boundary integral it turns out that the formulas for E^{z,w) 
(with the variables outside) and Hu{z,w) become exactly the same: 



Here one can easily identify the modulus and the angular part, for example 
for Ea{z,w): 



{z,w E C\Q). These formulas open up for geometrical interpretations. 

The exponential transform originally arose as a side product in the theory 
of hyponormal operators on a Hilbert space, see [5], [28], [30], and the inner 
product in the Hilbert space then automatically forces the transform to have 
certain positive definiteness properties, namely the two variable transforms 
1/Eq and Hq are positive definite (as kernels) in all C, and 1 — Eq is positive 
semidefinite outside fl. See [15] for direct proofs. For the four variable 
exponential transform we similarly have, for example, that 1/E is positive 
definite in the sense that 



for any finite tuples {zj,aj) e C^, e C, and with strict inequality unless 
all Xj are zero. 

Aside from operator theory, the ideas of the exponential transform appear 
implicitly as a tool in the theory of boundaries of analytic varieties, see [21], 






[3] (Chapter 19). 
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3 Harmonic, complex and exponential mo- 
ments 

The Cauchy and exponential transforms are generating functions for series 
of moments of a bounded domain Q G C Specifically we have the complex 
moments M^j = Mkj{Q), 



which in view of the Weierstrass approximation theorem completely deter- 
mine fl up to nuUsets, and the more restricted harmonic moments — 



The latter do not completely determine the domain, not even in the simply 
connected case (see [41]). However, for simply connected domains with ana- 
lytic boundaries they are at least sensible for local variations of the domain, 
ie., we have a local one-to-one correspondence 



A precise statement in this respect is the following: there exists a compact 
subset C r2 such that any Jordan domain D ^ K with Mk[D) = Mj^{Q) 
for all A; > necessarily agrees with Q. See [12], Corollary 3.10, for a proof. 

It is actually not a priori obvious that the moments {k > 0) are inde- 
pendent of each other (for example, the complex moments Mkj arc definitely 
not independent, given that they are moments of domains), but as will be 
discussed in Sections 7 and 8, it is indeed possible to define deformations of 
Q which change one of the without changing the other ones. See also [6] 
and [40] for further discussions and results. 

The generating properties for the harmonic and complex moments are 





(3.1) 



n ^ (Mo,Mi,M2,...). 



(3.2) 




(3.3) 



k=0 




k,j>0 
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Finally we introduce the exponential moments B^j — B^jip,) having 1 — 
E{z, w) as generating function, hence being defined by 

k,j>0 k,j>0 

The positive scmidefiniteness of 1 — E[z, w) implies that also the infinite ma- 
trix {Bkj) is positive semidefinite. The extended (four variable) exponential 
transform similarly generate some four index moments, see [32], [33], [36], 
[37], [28]. 

In the remaining part of this section we shall assume that is a simply 
connected domain with G and such that dfl is an analytic curve (without 
singularities). By writing (3.1) as a boundary integral. 



Mfc = ^ / z'' zdz, (3.5) 
Ida. 



27ri 



the harmonic moments become meaningful also for A; < 0, and these can be 
considered as functions of the moments for /c > 0, as far as local variations 
are concerned. Thus we may write 

M_k = M_fe(0) = M_fc(Mo, Ml, M2, . . . ) 

for A; > 0. 

In a series of papers by I. Krichever, A. Marshakov, M. Mineev-Weinstein, 
P. Wiegmann and A. Zabrodin, for example [53], [29], [27], [24], the so ex- 
tended moment sequence has been shown to enjoy remarkable integrability 
properties, for example (in present notation) 

This can be explained in terms of the presence of a prepotential J-'{fl) such 
that 

^ ,3.) 

In fact, for any sufficiently large > (it is enough that fl C D(0, R)), the 
energy functional 

J^{Q) = MQ) = \ ! [ log -^dm{z)dm{C) (3.8) 

^ Jn{o,R)\n Jo{o,R)\n l^~sl 
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serves as a prepotential. The exponential of can be identified with 

a certain r-function which appears as a partition function in mathematical 
models in statistical mechanics, see [23]. An indication of the proof of (3.7), 
and hence (3.6), will be given in Section 8. 

Remark 3.1. The dependence of M_fc and J^{Q) on the Mj is actually not 
analytic, so the above partial derivatives are really Wirtinger derivatives: 

d 1, d d . 



dM^ 2^dReMj dlmMj 

d 1, d d . 

+ 1- 



dMj 2'dReMj (9ImM/ 

On extending the power series (3.3) for the Cauchy transform to a full 
Laurent series by means of the negative moments one gets, at least formally, 
the Schwarz function of 

fc=— oo 

Even though the full series here need not converge anywhere it can be given 
the following precise meaning: the negative part of the series defines a germ of 
an analytic function at the origin and the positive part a germ of an analytic 
function at infinity. When is analytic the domains of analyticity of these 
two analytic functions overlap, and the overlap region contains dfl. Thus 
S{z), as written above, makes sense and is analytic in a neighborhood of dfl. 
Moreover, it satisfies 

S{z) = z for 2 e dn. (3.10) 

This can be seen by realizing that the positive part of the power series (3.9) 
defines the Cauchy transform Cfi{z), which can be written, for large z, 



Similarly, the negative part of the series agrees with the power series expan- 
sion of a boundary integral which can be interpreted as the Cauchy transform 
of the exterior domain: 

^ M, If (dC 
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for z small. Now (3.10) follows from a well-known jump formula for Cauchy 
integrals. 

4 Finitely determined domains 

In order to exhibit transparent links between some of the previously defined 
objects we shall now work in a basically algebraic framework. This will mean 
that we shall work with classes of domains involving only finitely many pa- 
rameters, specifically finitely determined domains in the terminology of [34], 
in other words quadrature domains [1], [42], [44], [17] or algebraic domains 
[48]. In this section our domains need not be simply connected, they will 
from outset be just bounded domains in the complex plane (bounded open 
set would work equally well). 

There are many definitions or equivalent characterizations of quadrature 
domains, the most straightforward in the context of the present paper per- 
haps being the following: a bounded domain fl in the complex plane is a 
(classical) quadrature domain (or algebraic domain), if the exterior Cauchy 
transform is a rational function, i.e., if there exists a rational function R{z) 
such that 

Ca{z) = R{z) for 8l\z&C\Vt. (4.1) 

Below is a list of equivalent requirements on a bounded domain C C. 
Strictly speaking, the last three items, c)-e), are insensitive for changes of 
by nuUsets, but one may achieve equivalence in the pointwise sense by 
requiring that the domain Q considered is complete with respect to area 
measure m, i.e., that all points a e C such that m(D(a, s) \ Q) = for some 
£ > have been adjoined to Q,. 

a) There exist finitely many points ajt e Q and coefficients c^j e C such 
that 

/ /i(im = J] J]c,.,/i(^)(a,) (4.2) 

•'^ fe=l j=0 

for every function h which is integrable and analytic in Q,. This is the 
original definition of quadrature domain, used in [1] for example. 

b) There exists a meromorphic function S{z) in Q, extending continuously 
to do. with 

S{z) ^z for ^ e dVL. (4.3) 
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This S{z) will be the Schwarz function of dQ [7], [44]. 



c) The exponential transform En{z,w) is, for z,w large, a rational func- 
tion of the form 

E„i.M-^^, (4.4) 
where P and Q are polynomials in one and two variables, respectively. 

d) Q is determined by a finite sequence of complex moments (Mjtj)o<fej<Ar 
(see Example 4.1 below for a clarification of the meaning of this). 

e) For some positive integer N there holds 

det(Sfej)o<jkj<iv — 0. 



Basic references for c)-e) are [34], [35], [36], [13]. 

When the above conditions hold then the minimum possible number N 
in d) and e) , the degree of P in c) and the number of poles (counting multi- 
plicities) of S{z) in Q, all coincide with the order of the quadrature domain, 
i.e., the number = XlfeLi^^ (4-2). For Q, see more precisely below. 

If fl is simply connected, the above conditions a)-e) are also equivalent 
to that any conformal map / : D — > is a rational function. If Q is multiply 
connected, then it is necessarily finitely connected (in fact, the boundary will, 
up to finitely many points, be exactly the algebraic curve Q{z, z) = with Q 
as in c) above), and there is a beautiful extension, due to D. Yakubovich, of 
the aforementioned statement: dfl is traced out (on 9D) by the eigenvalues 
of a rational normal matrix function F : D — )■ C^^^ which is holomorphic 
in D; see [54] for details. 

The positive definiteness properties of the exponential transform (see end 
of Section 2) imply that when i7 is a quadrature domain of order N then 
Q{z,w) admits the following representation [14]: 

Q{z, w) = P{z)PH - Yl Pk{z)P;M- 

k=0 

Here each Pk{z) is a polynomial of degree k exactly. Recall that P{z) = 
Pn{z) has degree A^, and one usually normalizes it to be monic (thereby also 
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making Q{z,w) uniquely determined). Thus the rational form (4.4) of the 
exponential transform can be expressed as 



N-l 



Pk{z)Pk{w) 
P(z)i>) 



With the same polynomials, the rational form (4.1) of the Cauchy transform 
is precisely 

Cn{z) - -p^. 

Note however that the double Cauchy transform Cn{z, w) will not be rational, 
only its exponential will be. 

Example 4.1. The following example should clarify the meaning of a domain 
Q, being finitely determined (determined by finitely many moments Mkj), as 
in the equivalent condition d) above. Suppose that, by some measurements, 
we happen to know the moments Mqo, Miq, Mqi, Mu, but we do not know 
from which domain they come. Generally speaking, there will be infinitely 
many domains whose moment sequence starts out with just these four mo- 
ments, unless of course the data are contradictory so that no domain has 
them. However, for exceptional choices of the data it happens that these 
four moments (or, in general, finitely many moments) determine the domain 
uniquely. This happens exactly for quadrature domains. 
Assume for example that 

Moo = Mio = Moi = 4, Mu = 12, (4.5) 

and let us try to compute the exponential moments Bkj as far as possible. 
We have 

1 - exp [ - V ^ifj' ] = l-exp[-^-^-^ - - ...] 

k,j>0 

4 4 4 12 1.4 ,2 1, 4 ,3 

= — H r H — H — H ( - ■■■) ( r -...)-... 

zw z'^w zw'^ z'^w'^ 2! zw 3! zw 

A A A A 

— H — — H — H — — - + higher order terms. 

zw z^w zw' Z'^UU^ 

Thus we see that 

Boo — Bio — Boi — Bii — A. 
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Clearly the determinant det(5jfcj)o<A;,j<Af vanishes for N — 1. The theory 
(equivalent condition e) above) then tells that the moments come from a 
quadrature domain of order one, i.e., a disk. To determine which disk one 
uses (3.4) together with (4.4). Arguing in general one first realizes that 

P{z)-PH E ^k^^i = Piz)PH - Qi^, w), (4.6) 

k,j>0 

in particular the left member is a polynomial. Setting P{z) = X^^q '^j^'' 
chosing ctTv = 1 (normalization) , the vanishing of the coefficient of 1/zw in 
the left member gives 

N 
k,j=0 

This, by the way, explains the condition det(i?fcj)o<A;j<iv = 0. It also gives 
the equation 

TV 

Y,Bkjaj = {k = 0,...,N) 

j=0 

for the coefficients of P{z). In our particular example this equation gives 
ao — — CKi = —1, hence 

P{z) ^z-1. 
Next, Q{z,w) is easily obtained from (4.6): 

Q{z, w) = P(.)PH - Piz)PH E ^^^1 

k,j>Q 

4 4 4 4 

= (ZW - Z - W + 1) - {ZW - Z - W + 1)(— + -ir-r H 3T + + . . . ) 

zw z^w zw^ z^w^ 
— ZW — z — w + 1 — 4 + terms containing negative powers of 2; or iD 

— ZW — z — id — 3 — {z — l){w — 1) — 4. 

Thus we have identified Q as the disk 

n = {zeC: |;2-1|^<4}. 

No other domain (or open set) has the first four moments given by (4.5). 
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5 The meromorphic resultant and the ehmi- 
nation function 



In this section we review the definitions of the meromorphic resuhant and 
the ehmination function, as introduced in [19], referring to that paper for 
any details. If / is a meromorphic function on any compact Riemann sur- 
face M we denote by (/) its divisor of zeros and poles, symbohcally (/) = 
/~^(0) — /~^(oo) (written in additive form). If D is any divisor and g is any 
function we denote by {D, g) the additive action of D on g, and by g{D) the 
multiplicative action of on D. See the example below for clarification. 
Now the meromorphic resultant between / and g is, by definition. 



whenever this makes sense (the resultant is undefined if expressions like O-oo, 
^, ^ appear when the expression in the middle member is evaluated). In 
the last expression, log (7 refers to arbitrarily chosen local branches of the 
logarithm. When defined, TZ{f,g) is either a complex number or 00. As a 
consequence of the Weil reciprocity theorem [51] the resultant is symmetric: 



Example 5.1. To illustrate and explain the above definitions we spell out 
what they look like for the sample divisor 



Here a,b,c G M (any compact Riemann surface) and g is a function on M 
(typically a meromorphic function, but this is not absolutely necessary as for 
the definitions). It is possible to think of D as a 0-chain and of as a 0-form. 
Then the above definitions amount to 



((/),log9> 



n{f,g)^n{g,f). 



D = 1- (a) + l- (6) - 2- (c). 




(5.1) 




(5.2) 



We can associate to D the 2-form current 



Sodx Ady — 6adx A dy + 6ijdx Ady — 26cdx A dy, 



(5.3) 
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where Sadx A dy denotes the Dirac measure (point mass) at a considered as a 
2-form (so that it acts on functions by integration over M). Then (5.1) can 
be augmented to 

{D,g) ^ / g ^ / gSodx A dy. 
Jd Jm 

If / is a meromorphic function with divisor (/) = D, then the above g{D) 
(in (5.2)) is also the resultant Tl{f,g) — g{{f))- For example, if M = P and 

{z-a){z-b) 
^^'^^ {z-cf ' 

then 

It is also possible to let functions of two (or more) variables act on divisors. 
For example, if g{z,w) is such a function and Di, D2 are two divisors, say 
= 1 • (a) + 1 • (6) - 2 • (c), ^2 = 3 • (p) - 3 ■ (g), then (by definition) 

g{Di,D2) =exp{Di(^D2,logg) = — .3 , r^. 5.4 

g{a,qrg{b,qrg{c,pr 

Next the elimination function between two meromorphic functions, / and 
gf, on M is defined by 

^f,gi^,w) ^ ^(/ - z,g-w), 

where z,w & C are parameters. It is always a rational function in z and w, 
more precisely of the form 

Q{z,w) 

^^-^'^""^ = Pizjm^y 

where Q, P and R are polynomials, and it embodies the necessary (since M 
is compact) polynomial relationship between / and g: 

SfMiO,9iO) = o (Cgm). 

There is also an extended elimination function, defined by 
Sf,g{z, w; a, h) = ^( J^^' 
15 



To relate the elimination function to the exponential transform one need 
integral formulas for the elimination function. For this purpose we shall 
make some computations within the framework of currents (distributional 
differential forms) on M. As a first issue, if / is a meromorphic function in 
M we need to make its logarithm a single- valued function almost everywhere 
on M. This is done as usual by making suitable 'cuts' on M, and by choosing 
a branch of the logarithm, call it Log /, on the remaining set M \ (cuts). Its 
differential in the sense of currents will then be of the form 

df 

dLogf = -J + distributional contributions on the cuts. (5.6) 

If there is a cut along the x-axis, for example, the distributional contribution 
will be dy times a measure on the cut, that measure being arc length measure 
times the size of the jump (an integer multiple of 27ri) over the cut, between 
the two branches of the logarithm. 
The next issue is the identity 

I ,,df, - , , 
— ci(y) = 5^f)dxAdy. 

See [19] for the simple proof. 

Now, proceeding formally we transform the resultant to an integral in the 
following series of steps: 

nL g) = 9{{f)) = e«^)'L°ss> ^ exp / Logy 

= exp[/ Log gd(f)dx A dy] = exp[— Loggd{—)] 

= exp — / — AdLogg\. 
27ri Jm f 

It should be remarked that, in the final integral, the only contributions come 
from the jumps of Log g (the last term in (5.6) when stated for Log g) because 
outside this set of discontinuities the integrand contains d( A d( = as a 
factor. Similarly one gets, for the elimination function, 

^fM, ^) = exp [— / A d Log {g - w)] . (5.7) 
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Now we shall apply (5.7) in the case that M is the Schottky double of 
a plane domain Q, and connect it to the exponential transform of Q. So 
let f2 be a finitely connected plane domain with analytic boundary or, more 
generally, a bordered Riemann surface and let 

M = n = nudnun 

be the Schottky double of Q, i.e., the compact Riemann surface obtained 
by completing fl with a backside with the opposite conformal structure, the 
two surfaces glued together along dQ (see [8] , for example) . On Q there is a 
natural anticonformal involution J : Q — > Q exchanging corresponding points 
on Q and Q and having dQ as fixed points. 

Let / and g be two meromorphic functions on Q. Then 

r-W^), 9*^(9^) 

are also meromorphic on Q. 

Theorem 5.1. With Q, Q, f , g as above, assume in addition that f has no 
poles in flu dfl and that g has no poles in QU dQ. Then, for large z, w, 

df 



/ 



z g* — w 



In particular. 



df 



A = 



df 



f-z f-w' 

Proof. For the divisors oi f — z and g — w we have, if z, w are large enough, 
supp(/ — z) <Z fl, supp(5f — w) G fl. Moreover, log(gf — w) has a single- valued 
branch. Log (g — w), in Q. Using that g = g* on dQ we therefore get, starting 
with (5.7), 

df 



^fA^^'^) = exp 



= exp 



= exp 



= exp 



exp 



1 

27ri 
1 

27ri 



f-z 

df 
f-z 



1 

27ri 



Ac? Log {g-w)] 
AdLog (g-w)] 
A Log {g - W)] 



df 



an f ~ z 

^ / ALog(^-w)] 
27n Jan f-z 



df 



A 



dg^ 



27ri / - ^ g 



w 



]■ 
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as claimed. 



□ 



We next specialize to the case that Q is a quadrature domain. Let S{() 
be the Schwarz function of Q. Then the relation (4.3) can be interpreted as 
saying that the pair of functions S{() and ( on fl combines into a meromor- 
phic function on the Schottky double fl of fl, namely that function g which 
equals S{() on Q, and equals ( on Q. 

The function f = g* = g o J is then represented by the opposite pair: ( 
on Q, S{() on Q. Therefore, on Q we have /(C) = g*{C) — C) Theorem 5.1 
immediately gives 

Corollary 5.1. For any quadrature domain fl, 

En{z,w) ^Sfj*{z,w) {\z\,\w\:^l), (5.8) 
where f , f* are the two meromorphic functions on Vl given on Vl by f{() = (, 

An alternative way of conceiving the corollary is to think of / being 
defined on an independent surface W, so that / : W — > Q is a conformal map. 
Then Q is a quadrature domain if and only if / extends to a meromorphic 
function of the Schottky double W, and the assertion of the corollary then 
is that the exponential transform of Q is given by (5.8), with the elimination 
function in the right member now taken in W. 

If fl is simply connected we may choose W = 3, so that W = F with 
involution J : ( Then / : D — >■ is a rational function when 

(and only when) Q is a quadrature domain, hence we conclude that Eci{z, w) 
in this case is the elimination function for two rational functions, /(C) and 

r(C) = 7(i7c). 

Finally we want to illustrate the effectiveness of the notation (5.4) by 
giving a formula for how the exponential transforms under rational conformal 
maps. 

Theorem 5.2. Let fl be bounded domain in the complex plane and let F 
be a rational function which is bounded and one-to-one on fl. Then for 
z,w & C \ F{fl) we have 

Ep^n^iz, w) = En{{F - z), {F - w)). (5.9) 
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Proof. We have 

i:. / N r 1 /" d(Ad( ^ 

F(n) 

F(a) 

= exp y" log(F(C) - ^) A rflog(F(C)-«;)j . 
n 

Let Dz — {F — z) denote the divisor of F{() — z, and consider it as a 
function : P — )■ Z, namely that function which at each point evaluates 
the order of the divisor (hence is zero at all but finitely many points). Then 

rflog(F(C) -z) = Y,D,{a)d\og{C - a), 

d log(F(C) - w) = Y.D^ (/3)c^ log(C - ;9) . 

With these observations we can continue the above series of equalities: 

^^(n) {z,w)^exp[^J d log(F(C) - ^) A (flog(F(C) - w)\ 

n 

a,/9eP 

which is (5.9). 

□ 
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6 Potential theoretic remarks 

Recall the definition (2.1) of logarithmic potential of a signed measure /j, with 
compact support in C. If /i has not zero net mass, then the behavior of the 
potential at infinity is such that it automatically puts an extra point mass at 
infinity, to the effect that the potential globally (on P) becomes the potential 
of a measure (again denoted fi) having vanishing net mass. If /j, is given by 
the divisor of a meromorphic (rational) function /, i.e., if 



The point we wish to make here is that the resultant TZ{f, g) between 
two rational functions, / and g, has a corresponding interpretation, namely 
in terms of the mutual energy. In general, the mutual energy between two 
signed measure fi and v of zero net mass (on P) is given by 



Some care is needed here since the above is not always well-defined (for 
example if ^ and v have common point masses). But this is just the same 
situation as for the resultant. Now, with / related as above to /i, and similarly 
g related to v hj dv = 5{^g)dx A d,y, then we have the relationship 



dji — S(^f)dx A dy 



in the notation (5.3), then 



U^^-log\f\. 





log I n{f, g) 



The proof is just a straight-forward computation: 



n{f,g)\^ = exp [((/), log^) + ((/),log^)] 




(/) 
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Similarly, one can relate the discriminant of / with a renormalized self- 
energy of fi. We refer to [19] for some details concerning this. 

All of the above generalizes well from P to an arbitrary compact Riemann 
surface M. The potential will then be less explicit, but it follows from clas- 
sical theory (for example Hodge theory) that given any signed measure /i on 
M with Jj^ d/i — there is potential C/^, uniquely defined up to an additive 
constant, such that, considering as a 2-form current, 

-d *dUf, = 271 fi. (6.1) 

Choosing in particular dfi = 6adx Ady — 6adx A dy for two points a,b E M 
gives a fundamental potential that we shall denote V{z) = V{z,w; a,b). 
Besides depending on a and b it depends on a parameter w, for normalization. 
The potential is characterized by the singularity structure 

V{z) — V{z, w; a,b) = — log |^ — a| -|- log |^ — 6| -|- harmonic, (6.2) 

together with the normalization V{w, w; a,b) — of the additive level. It has 
a few symmetries, 

V{z, w; a, b) = V{a, b; z, w) = —V{z, w; b, a), (6.3) 

and the transitivity property 

V{z, w; a, b) + V{z, w; b, c) = V{z, w; a, c). (6.4) 

If the Riemann surface happens to be symmetric (like the Schottky dou- 
ble), with involution J, then we have in addition 

V(z, w; a, b) = V(J(z), J{w); J{a), J(b)). (6.5) 

As an application, the ordinary Green's function g{z, Q for a plane domain 
(or bordered Riemann surface) Q, can be obtained from the potential V for 
M = Qby 

g(z, C) = lv(z, J(z); C, J(C)) = V(z, w; C, J(C)), (6-6) 

where w is an arbitrary point on dfl and where the second equality follows 
from (6.3), (6.4), (6.5). Cf. [9], p. 125f. See also [16] for some more details. 
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The potential V appears implicitly or explicitly in classical texts, e.g., 
[52], [43]. Here we discuss it briefly because it connects to the exponential 
transform. To start, it has a certain self-reproducing property, expressed in 
the identity 

V{z, w;a,b)^ — [ dV{-, c; z, w) A* dV{-, c; a, b), (6.7) 
27r Jm 

valid for arbitrary non-coinciding points z, w, a,b,c & M. Alternatively, 

V{z,w;a,b) = - [ dV{-,c;z,w) AdV{-,c;a,b), (6.8) 

where dV = ^dz, dV — ^dz. The proof of (6.7) is straight-forward: the 
right member becomes, after partial integration, an application of (6.1) and 
by using the symmetries (6.3) and (6.4), 



-— / V{-,c]z,'w) /\d* dV{-,c]a,b) 
27r Jm 



— V{a, c; z, w) — V{b, c; z, w) — V{z, w; a, b), 

which is the left member. 

In the case of the Riemann sphere, V relates to the cross-ratio (2.7 by 

Viz, w;a,b) = -\og\{z:w:a:b)\^-log\ \' ~ ~ ^\ \ . (6.9) 

^ [z — b){w — a)^ 

Therefore the extended exponential transform for any domain Q C P can be 
conveniently expressed in terms of V for P, for example 

Eq{z, w; a, b) = exp [ — / dV{-, c; z, a) A dV{-, c; w, 6)1 , 
ITT Jn 

and for the modulus, 

\EQ{z,w;a,b)\ = exp \ / dV{-,c; z,a) A *dV{-,c;w,b)]. 

Jn 

In the case of a compact Riemann surface, Jj^ dV{-, c; z, a) A dV{-, c; w, b) 
is purely imaginary, by (6.8). In particular we get, for Q = M = P, and by 
using (6.8) and (6.9), that 

Efiz, w; a, b) = exp / dV{-, c; z, a) A dV{-, c; w, b)] 

ITT Jp 
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= exp [ - 2V{z,a;w,b)] ^ \{z : a : w : b)f . 

Thus formula (2.6), which was asserted and proved in another way in Sec- 
tion 2, has been proved again. 



7 Moment coordinates 

In this section we restrict again to the case that O C C is a simply connected 
domain with analytic boundary, and with G ^2. In order to keep the 
presentation reasonably rigorous we shall, much of the time, actually restrict 
even more, namely by considering only those domains which are conformal 
images of the unit disk under polynomial conformal maps. 



with Oq > (normalization) and N fixed. This will be a manifold, call it 
M.2N+1, of real dimension 2N + 1, in fact it is an open subset of M x C'^ = 
R2^+^ with ao e K and e C (1 < J < N) as coordinates. 

Conformal images of the unit disk under polynomial conformal maps are 
exactly the simply connected quadrature domains with the origin as the sole 
quadrature node, i.e., with m = 1, ai = in the notation of (4.2), and the 
order of the quadrature domain then equals the degree of the polynomial 
{N + 1 in the above notation). The class of domains we are considering 
here is slightly smaller than this class of quadrature domains because the 
assumption of analytic boundary require our polynomials / to be univalent 
in a full neighborhood of the closed unit disk. But it has the advantage 
that M.2N+1 will be an ordinary open manifold (without boundary). For 
algebraic purposes one could equally well work with the just locally univalent 
polynomials, i.e., those polynomials / in (7.1) for which /' 7^ on D. 

Remark 7.1. The restriction to polynomial / does not change anything in 
principle, all formulas will look basically the same as in the general case. In 
the other extreme, without changing much in practice one may pass from 
the infinite dimensional case to the formal level (with no bothering about 
convergence) of transfinite functions, see [18]. 

As real coordinates on M.2N+1 we may of course use the real and imag- 
inary parts of the coefficients appearing in (7.1), namely Qq, ReUj, Imuj 



N 




(7.1) 



k=0 
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(1 < j < N), but much of what we are going to discuss concern what hap- 
pens when we switch to Mq, KeMj, ImMj (1 < j < A^) as local coordinates. 

The map (or change of coordinates) (oq, ai, 02, . . . ) ^ {Mq, Mi, M2, . . . ) 
is explicitly given by Richardson's formula [39], 



The summation here goes over all multi- indices (jq, • • • ,jk) > (0, • ... 0) for 
which jo + ■ ■ ■ + jk + k < N. It is convenient to set aj = for j > N, and 
then the upper bound is not needed. (And with Uj = for j < the lower 
bound is not needed either.) The formula (7.2) is proved by pulling (3.5) 
back to the unit circle by the map (7.1) and then using the residue theorem 
(after having replaced / by /*). 

It is clear from (7.2) that M^r+i = Mjv+2 = • • • = when cat+i = 
ciN+2 — • • • = 0, and it follows from basic facts about quadrature domains 
that the converse statement holds as well (when / is univalent in a full 
neighborhood of the closed unit disk). It is virtually impossible to invert 
(7.2) in any explicit way. Even for = 1 one would need to solve a third 
degree algebraic equation in order to express Oq, ai in terms of Mq, Mi. And 
the map (ao, oi, 02, . . . ) ^ (Mq, Mi, M2, . . . ) is (generally speaking) not one- 
to-one, as the example [41] (and several similar examples, cf. [55]) shows. 
The author does not know, however, of any such example in the present 
setting of polynomial conformal maps. C. Ullemar [47] gives an example of 
two polynomials of degree three, one univalent and one only locally univalent, 
which give rise to the same moments. 

Remark 7.2. Even when Mjv+i = Mi\i^2 = ■ ■ ■ = 0, the negative moments 
M_i, M_2, . . .generally make up a full infinite sequence (of nonzero num- 
bers). This is related to the fact that when Q is a quadrature domain, then 
its complement is almost never a quadrature domain (one can easily make 
sense to the notion of an unbounded quadrature domain, or a quadrature 
domain in P). 

A deformation of a domain Q as above corresponds to a smooth curve 
t fl{t) with ^2(0) = fl (equivalently, t ^ /(C) ^)) iii ■M.2N+1, and its 
velocity at i = is a vector in the real tangent space of M2N+1 at the point 
in question. This tangent space has, in terms of the coordinates introduced 




(7.2) 
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above, two natural bases, namely 

d d d 

duQ ' 9 Re ' d Im aj 

d d d 

aMo' dReMj' dlmMj 

(1 < j < A^), respectively. 

It is natural to consider also the corresponding complexified tangent 
space, obtained by allowing complex coefficients in front of the above ba- 
sis vectors. For this tangent space we also have the bases 

d d d 

d d d 

9Mo' dWj' Wj' 

Similar considerations apply for the cotangent spaces, i.e., we have the two 
natural bases, dao, daj, ddj and dMo, dMj, dMj, etc. It should be noted, 
however, that complex tangent vectors such as (j > 1) do not really 
correspond to velocity vectors for deformations of domains, only vectors in 
the real tangent space do. 

The dependence of / on the coefficients ao,...,aN is certainly analytic 
by (7.1), but as can be understood from the appearance of conjugations in 
(7.2), the dependence on Mq, . . . , Mjv is no longer analytic. Therefore we 
prefer to write this dependence as 

/(C) = /(C; Mjv, . . . , Ml, Mo, Ml, ... , M^), (7.3) 

or just /(C; M), with M shorthand for the hst of moments: 

M = (M^, . . . , Ml, Mo, Ml, ... , M^). (7.4) 

Despite the implicit nature of the dependence of / on the moments, sev- 
eral interesting general statements can be made, and we want to highlight 
a couple of them. The first is due to O. Kuznetsova and V. Tkachcv [25], 
[46], who were able to compute the Jacobian determinant for the transition 
between the above sets of coordinates, thereby confirming a conjecture of 
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C. UUemar [47]. The beautiful result can be expressed in terms of the ratio 
between the volume forms as follows: 



dMN A . . . dMi A dMo A dM^ A • • • A dM^ = (7.5) 
=2a^'+2^+^ 7^(/', /'*) -ddN ^...ddl^daQ^dal^■■■ ^ dttN- (7.6) 

Here TZ{f', /'*) is the meromorphic resultant between /' and /'*, as discussed 
in Section 5. When / is univalent (or locally univalent), then TZ{f', /'*) 7^ 0, 
hence the coordinate transition is locally one-to-one. 

The second general statement is actually a series of results obtained by 
M. Mineev-Weinstein, P. Wiegmann, A. Zabrodin, with later contributations 
also from other authors, [53], [29], [27], [24]. We shall below, and in the next 
two sections, discuss some of these results, namely the string equation (7.8), 
the prepotential (3.8) and integrability of moments (3.6) (we will provide a 
proof of (3.7)) and, finally, Hamiltonian formulations of the evolution equa- 
tions for the mapping functions (9.8). 

The results arc partly formulated in terms of a Poisson bracket, which for 
any two functions / = f{(; M), g = g{(; M), analytic for ^ in a neighborhood 
of the unit circle and also depending on the moments M, is defined by 

|/,^|-C^^-C^^- (7.7) 

As a first choice here we take / to be the conformal map (7.1) and g — f*, 
the reflection of / in the unit circle: f*{(;M) = f{l/(;M). Since we are 
assuming that Q has analytic boundary, / and /* are both analytic in some 
neighborhood of 9D, hence so is {/,/*}. Now it turns out that it is much 
better than so, in fact {/, /*} is analytic in the whole Riemann sphere, hence 
is constant. This remarkable fact is sometimes referred to as the string 
equation, which more precisely reads 

{/,/*} = !• (7.8) 

As for the proof (which follows [53] etc.) one first notices that f* = Sof 
{S is the Schwarz function), since this holds on dB>. Since also S — S{z; M) 
depends on the moments this spells out to 

r(C;M)=5(/(C;M);M). 
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Using hence the chain rule when computing one arrives at 




of). 



Now we notice from (3.9) that 




— h positive powers of z, 



since the coefficients Mi,M2, ... of the remaining negative powers of z in 
(3.9) are independent of Mq. From this one sees that {/, /*} is holomorphic 
in D. Similarly, it is holomorphic in P\ D, hence it must in fact be constant, 
and one easily checks that this constant is one, proving (7.8). 

The string equation (7.8) may appear to be quite remarkable identity for 
conformal maps, but actually it can be identified with an equation which 
was discovered much earlier, namely the Poluharinova-Galin equation [10], 
[31], [50], [20], [49] for the evolution of a Hele-Shaw blob with a source at the 
origin (Laplacian growth), this taken in combination with the Richardson 
moment conservation law [39] for such a flow. These considerations will be 
elaborated in the next section. 



8 Domain deformations driven by harmonic 
gradients 



Continuing on the theme of the previous section, we now discuss general 
domain variations from the point of view of Laplacian growth. We first con- 
sider an, in principle arbitrary, deformation t i-> Vt{t) of a simply connected 
domain Q = Q(0). Let Vn denote the speed of V{t) — dQ{t) measured in the 
outward normal direction. We need to assume that all data are real analytic. 
This can for example be expressed by saying that the Schwarz function S{z, t) 
of T(t) is real analytic in the parameter t (besides being analytic in z, in a 
neighborhood of T{t)). Then, by the Cauchy-Kovalevskaya theorem, there 
exists, for each t, a harmonic function g{-,t) defined in some neighborhood 
of r{t) and satisfying 




on r(t) 
on r{t). 



(8.1) 
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This function g can actually also be defined directly. In fact, it is easy to see 
that (8.1) can be expressed as 




holding on r(t), and hence identically (both members being analytic in z). 
To be precise, the latter equation actually contains (8.1) only with the first 
condition weakened to g{-,t) being constant on r{t). Prom the above one 
gets the formula 



where ^o(^) is an arbitrary point on r{t). Certainly, there will be exist 
neighborhood of r(0) in which all the functions g{-,t), for \t\ sufficiently 
small, are defined. 

Near r(0) we have the usual cartesian coordinates {x, y) {z — x + iy), but 
from g we also get another pair, {t, 9), namely 



Here *g is the conjugate harmonic function of g. We shall also write d6 = 
— *dg, so this star is actually the Hodge star. Thus t = t{z) is the time when 
r{t) arrives at the point z, and ^ is a coordinate along each T{t). In order 
for t{z) to be well-defined, at least locally, we need to assume that T4 0- 

Fix a value of t and introduce auxiliary coordinates (n, s) in a neighbor- 
hood of r(t) such that, on r(t), ^ is the tangent vector and is the outward 
normal vector (thus s is arc length along T{t)). Then we have (on T{t)) 




t = t{z) determined by z e ^t{z), 
0^-*g{z,t{z)). 



*dg — 




and so 



* dg{z,t) = Vnds. 



On the other hand (still on r{t)), 



— -0 — 
ds ^ dn Vn 
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so that 



dt = —dn. 



n 



It foUows that dt f\dO = dn f\ ds on V{t). But since ^) is just a rotation 
of we also have dn A ds — dx A dy. This gives dt A dO — dx A dy 

on r(i), which does not involve the auxihary coordinates {n,s). Hence the 
latter equation actually holds in a full neighborhood of T{t). Equivalently, 



We summarize. 

Proposition 8.1. Any real analytic deformation t ^ T{t) can be described, 
locally in regions where 14 7^ and in terms of local coordinates {t, 6) intro- 
duced as above, by the "string equation" 



Compare discussions in, for example, [22]. The reason for calling (8.3) a 
string equation will become clear below. 

Now we change the point of view. In the above considerations, the normal 
velocity of the boundary was the given data (essentially arbitrary, except 
for regularity requirements) and the harmonic function g was introduced as a 
secondary object. Laplacian growth, on the other hand, starts with a rule for 
prescribing a harmonic function g = g^, associated to any given domain fl, 
defined in the domain (or at least near the boundary of it) and vanishing on 
the boundary. Then one asks the boundary to move with the normal velocity 



Thus g now works as the generating function which forces the motion of 
dQ, and setting again 9 = — * g, (8.3) remains valid in regions where Vn 7^ 
0. Typically, g is specified by suitable source terms in Q or by boundary 
conditions on some fixed boundaries. 

The classical case is that e Q and that g = ga is the Green's function 
of Q with pole at 2; = 0: 



d{x, y) 
dit, 9) 



= 1. 




dt Ad9 = dx A dy. 





g{z, t) = — log 1^1 + harmonic 




(8.4) 
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When pulled back to the unit disk by / : D — )■ f2 (normalized conformal map) 
one then gets g(f((^),t) = — log \(\, and hence 6 = —*g{f{C),t) = arg(^. Thus 
9 will simply be the angular variable on the unit circle. In addition, one sees 
that d6 can be interpreted as harmonic measure (on 5D, and on dQ). This 
classical version of Laplacian growth is connected to the harmonic moments 
in that it is characterized by the moments M^, k > 1, being constants of 
motions ('integrals'), while Mq evolves linearly in time. Precisely, with g 
defined by (8.4) we will have (as functions of t) 

Mo — 2t + constant, 
Mk — constant [k > 1). 

These conservation laws were discovered by S. Richardson [39], and they are 
easy to verify directly. In fact, one shows that 



d 
~dt 



hdm = 27r/i(0) 

Q(t) 



for any function h which is harmonic in a region containing the closure of 
n(t), and then (8.5) is obtained by choosing h{z) = (k > 0). 

Now we can compare (7.8) with (8.2) or (8.3): since {/, /*} is an analytic 
function, the restriction of (7.8) to 83 contains all information. And on dB> 
we have, on setting ( = e'^ and writing /(e'^, t) = x + iy, 

{/,r} = 21m[^#^]=2-^^"'^) 



■dOdMo^ d{Mo,e)' 

Since 2^ = f by (8.5) it follows that (7.8) is the same as (8.2), or (8.3). 
This was the reason for caUing (8.3) a string equation. 

Renicirk 8.1. A common variant of the above is that Q is a domain in the 
Riemann sphere, oo e Jl, and that g is the Green's function with pole at 
infinity. 

Remark 8.2. Most of the above considerations also apply in higher dimen- 
sions, the main difference being that, in n real dimensions, dO — — * dg will 
be an (n — l)-form and that (8.3) will take the form dt Ad9 = dxi A - ■ ■ A dxn 
(the volume form). Another generalization, elliptic growth, is discussed in 
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Now we shall connect Laplacian growth evolutions to the derivatives -^-r, 



dM. 
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These evolutions will not be monotonic, so we must allow now Vn to 

change sign, and hence to be zero at some points of the boundaries. Let 
L denote a real-valued differential operator with constant coefficients acting 
on functions in C, for example L — Then, given f2(0) with analytic 

boundary, there exists locally an evolution such that 



~dt 



I hdm^27r{Lh){0) (8.6) 
Jn{t) 

for every function h which is harmonic in a region containing the closure 
of Q{t). This is a version of Laplacian growth, in fact it is the evolution is 
driven by the harmonic function 

giiz) = {Lag{z,a)}a=o, 

where g{z, a) — — log — a| + harmonic is the ordinary Green's function 
with pole at a, and the differential operator = L acts on that variable (the 
location of the pole). Thus, gi is harmonic in fl, has a multipole singularity 
at the origin and vanishes on dQ. 

It is actually not easy to prove that evolutions fl{t) driven by a given 
harmonic gradient really exist, but it can be done under the right assumptions 
(analyticity of the initial boundary is usually necessary, for example). An 
example of such a proof is given in [38]. Other approaches are presented in 
[45] and [26]. In the algebraic framework of the previous section, i.e., with 
fl G A^2Af+i, the full evolution will stay in M.2N+1 provided the order k of 
the differential operator L is at most N, and the existence proof is much 
easier, see [11]. In fact, the existence and uniqueness can even be read off 
from (7.5). 

For the rest of this section, and also for the next section, even if it does 
not change anything in practice it is useful to think of being in the finite 
dimensional setting of A42N+i, hence assuming that k < N. 

To prove the assertion that gL really achieves (8.6), recall that the evo- 
lution t ^ fl(t) is defined by the outward normal velocity of the boundary 
being Vn — — ^^gjf^ ■ Since —Ag{-, a) = 2Ti5a we have 

-Ag^L = -{AL„g'(-, a)}a=o = -{LaAg{-, a)}a=o = 27rL5o. 

With h harmonic this gives 

^ f U J f TA L J f u^9l 



I hdm— I Vnhds — — f h- 
dt Jn(^t) JdQ.(t) Jan{t) c/n 



ds 
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= - / hAgL dm = 27rL/i(0), 
Jn{t) 

proving (8.6). 

Useful particular choices of L above are, for A; > 1, 

Li 



We may allow h in (8.6) to be complex- valued, and choosing h{z) — and 
evaluating at 2; = one gets 

{Li(^'=)},=o = fc!, {L2(^')}.=o = i • fc!, (8.7) 

and {Li{z^)}z=Q = (i = 1,2) in all remaining cases. Similarly, with h[z) = 
z^ we have 

{i^i(^')}.=o = A;!, {L2(^'=)}.=o = -i • A;!, (8.8) 

and {Li{z^ )} z=o = in the remaining cases. 

In what follows, we will denote the time coordinates for the two evolutions 
corresponding to Li and L2 by different letters, namely ti and respectively. 
It is a consequence of (8.6) and (8.7), (8.8) that the evolution ti 1-^ ^{t\) 
generated by ql^ has the property that 



Similarly, for the evolution generated by ql^: 



In other words, the tangent vector giving the velocity for the evolution driven 
by Qli is 

d d d „ , , d 

dh '^dMk dMk 'dReMk 

and that for 51^2 is 

d , d d . , , d 
— = 2ik\{—— - — ^) = 2k\- 



dt2 '^dMk dMk' 'dlmMk 

Now we wish to extract jSr- and from the above relations. We 
then see that these partial derivatives will correspond to linear combinations 
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with complex (and non-real) coefficients of the two different evolutions, or 

tangent vectors, hence they will not in themselves correspond to evolutions 
of domains. Still they make sense, of course, as vectors in the complexified 
tangent space of A^27v+i- Precisely, we get 

d 1 . d . d , 



dMh AkVdti dt2' 
d I , d . d ^ 

= 7r7(3r + iir)- 



dMk AM^dti dt2' 
This is what we need in order to prove (3.7). Using (3.8) we get 

dF{n) d 1 



dMf. OMkTT Jo{0,R)\Q. Jo{o,R)\n \^~Q\ 



Air^kl'dti 



' ^ '^'ir^ I I log I ^ dm{z)dm{C,). 

"'^2 Jn{Q,R)\n JB{o,R)\n \^ ~ s\ 



When computing the for example, Q is to be replaced with the evolution 
fl{ti), and we evaluate the derivative at ti = {fi{0) = Q). There are two 
occurrences of Q in the expression for J-{fl), and ^-"(^2) symmetric in them, 
so it is enough to differentiate one of the occurrences and then multiply the 
result by two. Writing 2 log l-z — CI — log(C ~ ^)+ log(C ~ using (8.6), 



(3.5) this gives 

r\u=o f f log\z - C\dm{z)dm{C) 

h J'D(o,R)\n{ti) JB{o,R)\n{ti) 

A| ( f - f {logic -z)+\ogiC-z))dm{z))dmiC) 
"ii J»(Q,R) Jn(u) 



d 

dU 



I 

Jb( 



>,R)\n "f-i J»{o,R) Jn{ti) 

Jo{o,R)\n ^ 



- - 1)1 £ (^ + ^) CdC = -27r^^!(iM_, + iM_,). 
Similarly, the corresponding derivative with respect to t2 gives 



d I 



D(o,ii)\n(t 



/ log 

i) Jn{o,R)\n{ti) 



z — C\ dm{z)dm{() 
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Hence, in the combination ^ — i^, M_fe will cancel out, and we get 

dj^m _ 1^ 

as desired. 



9 Hamiltonians 

As a final topic, we wish to explain, in our notations and settings, the Hamil- 
tonian descriptions of the evolution of the mapping function / presented in 
[53] and related articles. Recall that we work within M.2N+I: and we use 
M as a shorthand notation for the moments, which serve as coordinates on 
M2N+1, see (7.4). 

Let W = W{z; M) be a primitive function of the Schwarz function S — 
S{z; M). By (3.9), the power series expansion of W will be 

W{z-M)^- J2 + Mo log z + C(M), (9.1) 
kei\{o} 

where C(M) is a constant. Decomposing W into real and imaginary parts, 
W — w + i * w, the real part can be made perfectly well-defined. In fact, 
S{z) — 2^, and since S{z) — z on dQ we see that the real- valued function 

1 1 

u(z) = -\z\'^ w(z) 

satisfies ^ — on dQ, hence u is constant on dQ. We shall fix the free 
additive constant in u, and hence that in w, by requiring that it = on dQ,. 
It follows that u then satisfies 

in a neighborhood of 9Q, 
on dn. 

Next we note from (9.1) that 




^ ifc>0 ^ " 
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Thus, as for the real part, = log \z\ + harmonic in Q. On dfl we have 
~ ~'^'§Mi ~ ^' because u vanishes to the second order on Hence it 
follows that is minus the Green's function g — go^oi^ with pole at the 
origin: 

Let G = gi + i * 51 be the analytic completion of g. This Green's function 
G — G{z; M) is (hke g) a function of 2; e Q and the moments. By (9.2), 

Recall that / = /(C; M) denotes the conformal map (7.1) from D to Vt. Define 
the (complex- valued) Hamiltonian function Hq = Hq{(; M) associated to Mq 

by 

i/o(C; M) = G(/(C; M); M) = - log C, (9.4) 

the last equahty because G o / is the Green's function of D. We see that Hq 
actually is independent of the moments. Therefore, as |^ = 0, — ~1 
by (9.4), we trivially arrive at the evolution equation 

W,Ho}, (9.5) 



9Mo 



valid for any smooth function ip = ip{(] M). 

For the higher order moments, there are corresponding identities, and 
they are more interesting and more selective: they hold only for (p — f. We 
first define, for /c > 1, 

H,{C, M) = - ^^^^^\ where z = f{C,M). (9.6) 

The above means that the derivative only acts on the which appears 
in W{z; M), not that in f{(; M). Always when we write z in place of /(C; M) 
it shall have this implication. Now, G{f{(; M); M) — — logC gives again, by 
differentiating with respect to ( and M^, 

^(.M)|(C;M) = -i, 



35 



f)C r) f r)C 

^(.;M)3^(C;M) + ^(.;M))=0. 

Thus, multiplying the latter equation with C,^ and using the first equation 
gives 

^(C; M) = C^(C; M)^(z; M). (9.7) 

Here we wish to remove G the from the right member, in favor of H]^. 
Prom (9.6) we get 

<t(C;M)^-^(.M).c|(CM) 

and, using also (9.3), 

The last term coincides with the right member in (9.7). Thus substituting 
into (9.7) we get 

= (9.8) 

which is the evolution equation for / we wanted to arrive at. 

As for the classical Laplacian growth evolution, even though the situation 
is essentially trivial in the moment coordinates (in view of the explicit solu- 
tion (8.5)), it may be nice to put it all in a traditional Hamiltonian framework 
(with real- valued Hamiltonian function). One possibility then is to choose 
the phase space to be (9D x A^2Af+i C x M^^+^, with real coordinates 
(6*, Mo, Re Ml, Im Ml ... , Re Mat, ImMjv), symplectic form 

CO =^d9 A dMo + dReMi A dImMi H h dReMjv A dlmMjv 

=^ [ide A dMo + dMi A dMi H h dM^v A dM^] 

and Hamiltonian function 

ii'(^,Mo,ReMi,...,ImMjv) = 6. 
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The Hamilton equations are, generally speaking (see [4]), 

-dH = co{C,-), (9.9) 

where C = ^ is the velocity vector for the evolution, a vector in the tangent 
space of the phase space. In our case we have (using dot for time derivative) 

■d-d ■ d ■ d 

i = 6'— + Mo—— + Re Ml H hlmM^r- 



86 ^OMq aReMi aimMjv' 

giving in (9.9) 

^ = 0, Mo = 2, 

ReMj = ImM,- = (1 < j < AT), 

as expected (cf. (8.5)). 

Note that the first term in tj, with dM^ = 2dt, can be identified with 
—dx A dy, in the notation of Proposition 8.1. 
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